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Systems of exchange-coupled spins are commonly used to model ferromagnets. The quantum 
correlations in such magnets are studied using tools from quantum information theory. Isotropic 
ferromagnets are shown to possess a universal low-temperature density matrix which precludes en- 
tanglement between spins, and the mechanism of entanglement cancellation is investigated, revealing 
a core of states resistant to pairwise entanglement cancellation. Numerical studies of one-, two-, and 
three-dimensional lattices as well as irregular geometries showed no entanglement in ferromagnets 
at any temperature or magnetic field strength. 
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PACS numbers: 75.10. Jm,03.67.-a,03.67.Mn 



Quantum correlations are responsible for much of the 
challenge in understanding interacting many-body quan- 
tum systems, and it is therefore of fundamental impor- 
tance to have quantitative knowledge of these correla- 
tions. Progress in quantum information theory has led 
to the development of new measures of the inseparabil- 
ity of a quantum state, and in the last few years these 
measures have been used to assess the quantum correla- 
tions in diverse physical systems. Concurrence is an 
especially useful metric for such studies because it can 
be applied to mixed as well as pure states. It therefore 
can be used to quantitate the thermal entanglement in 
a system at nonzero temperature. It can also be applied 
to evaluate the inseparability of an equal incoherent mix- 
ture of degenerate energy eigenstates. However, concur- 
rence is defined only for a pair of qubits. Since a qubit is 
formally equivalent to a spin-1/2 particle when only the 
spin degree of freedom of the latter is considered, this has 
led to several analyses of the thermal entanglement be- 
tween a pair of interacting spin-1/2 particles. Two-spin 
systems studied in this way include the XXX 0, XX Z 
[g, and XYZ |i Heisenberg models, the XX-model 0, 
the XF-model [Q], the Heisenberg-Dzyaloshinski-Moriya 
model , and the Ising model [J Inseparability at zero 
and finite temperture has also been studied in two-spin 
subsystems of one-dimensional spin chains, including the 
XXX-chain 0,0,0, the XF-chain the Ising 

chain 0, and the Majumdar- Ghosh chain as well as 
in the ground state of quasi-one-dimensional spin ladders 

This paper will examine quantum correlations in 
isotropic ferromagnets without restrictions on the geom- 
etry of the spins (such as one-dimensionality or period- 
icity) or on the range of the spin-spin interaction. A 
suitably general model of a magnet, which will be called 
the Heisenberg spin graph, consists of N spins coupled by 
exchange interactions with arbitrary range and strength. 
It is defined by the Hamiltonian 



This can be visualized as a graph with a spin-1/2 parti- 
cle at each vertex, where the edge between spins i and 
J, representing the exchange interaction • Sj between 
the two spins, is weighted by the coupling constant J^ . 
The summation runs over each pair of spins once. The 
operator Si — {Sf,Sf,S^) is the spin operator associ- 
ated with the spin at the ith vertex of the graph; its 
Cartesian components obey the usual angular momen- 
tum algebra. As a matter of notation, states in the 2^- 
dimensional Hilbert space H of the model are specified 
with respect to basis states of the form \i,j,...,k), where 
the listed spins i,j,...,k are aligned in the -|-z-direction 
('up') and the remaining spins are antiparallel ('down'). 
|0) denotes the state with all spins down. One might 
anticipate that the Heisenberg spin graph would be so 
general a model that nothing much could be said about 
it. It turns out that when the coupling constants Jy are 
nonpositive (that is, ferromagnetic), certain features of 
Heisenberg spin graphs are universal, in the sense that 
they do not depend on the exact values of the coupling 
constants or even on the number of spins. This univer- 
sality allows a complete analysis of low-temperature cor- 
relations in such ferromagnets. 

Because the Heisenberg spin graph Hamiltonian Hhsg 
commutes with the total 2:-spin operator 5^ = 
the two-qubit reduced density matrix for any two spins in 
an eigenstate of Hhsg always takes the restricted form 
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Hhsg = ^ JijSi • Sj 



(1) 



For such a density matrix, the concurrence is C(p) — 
2 max (0, 17I — \/ae). Concurrence ranges from zero, for a 
separable state, to one, for a maximally entangled state. 
The entanglement of formation, another important en- 
tanglement measure, can be calculated directly from the 
concurrence, and is monotonically related to it. The 
method of calculating the concurrence for more general 
density matrices can be found in Wootters . 

In Ref. the authors have reported the absence of 
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thermal and magnetic entanglement in the ferromagnetic 
Heisenberg spin chain, based on the analytical solution 
of the N — 2 case and numerical studies of < 10 
chains. The following argument shows that the absence 
of thermal entanglement, at least at low temperatures, 
is a general feature of exchange-coupled ferromagnets, 
irrespective of the range of the exchange interactions, 
the dimensionality of the solid, or the precise bound- 
ary conditions (free or periodic), even in the presence 
of inhomogeneous coupling strengths. Consider the fer- 
romagnetic Heisenberg spin graph Hamiltonian for N 
spins: Hhsg = Y^i<j Jij^i ■ Sj: where J^j < 0. Be- 
cause the exchange interaction can be rewritten in terms 
of a permutation operator and the identity operator as 
Si • = i (Vij — iX) , the N + 1 completely symmetric 
states 
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are eigenstates of Hhsg with the common eigenvalue 
i Si< j '^ij ■ These states span the ground subspace of 
the model; there are no 'accidental' degeneracies (see 
Appendix). At sufficiently low temperature, the thermal 
density matrix is an incoherent mixture of these N + 1 
states, the contribution from other states being negligi- 
ble. 

We choose any two spins in the solid and trace out the 
rest. Because of the linearity of the partial trace, we can 
perform the trace on each of the completely symmetric 
states separately. The total reduced density matrix p will 
then be 
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where p„, the reduced density matrix of the nth com- 
pletely symmetric state, has the form: 
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From combinatoric reasoning we have 
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Pn = 7n = Sn = 
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for the components of the density matrix p„. 

Carrying out the summation we find that the ground 
state density matrix is 
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Thus the low temperature density matrix for any two 
spins assumes a universal form, irrespective of the 
strength, geometry, or homogeneity of the interactions 
or the dimensionality of the solid. The concurrence of 
this universal density matrix is zero. This shows that, at 
least at low temperatures, no two qubits in an isotropic 
ferromagnet share any entanglement. The correlations 
present, such as {SfSJ) — 1/12, can be regarded as 
purely classical. 

The form of the reduced density matrix given in Eq. @ 
is of course strictly valid only in the limit of zero tem- 
perature. As the temperature T increases, the remaining 
eigenstates of Hhsg (those not listed in Eq. Q) begin to 
contribute to the reduced density matrix p describing the 
state of the two chosen spins. As there are a finite number 
of such states, it can be seen that the elements a, /?, . . ., 
of the reduced density matrix are continuous functions 
of T. Since at zero temperature I7I — ^/ae — —1/6, there 
must be some (possibly infinite) temperature interval, 
with T = as one endpoint and closed if the interval is 
finite, in which the concurrence is zero. Thus the onset of 
entanglement in ferromagnets, if it occurs at all, happens 
as a phase transition at nonzero temperature. 

Note that, individually, each p„ is an inseparable state, 
except for po ^-nd pn, with concurrence 



C{Pn) = 



N{N - 1) 
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(n{N - n) - ^n{n - - n){N -n - 1)) . 

It is the incoherent mixing of these inseparable states 
that washes out the quantum correlations. One might an- 
ticipate that this entanglement cancellation occurs pair- 
wise, between states related by flipping the orientation 



of each spin, so that p'„ = (p,, 



/2 would be sep- 



arable; only the entanglement of Pn/2 when N is even 
could not be eliminated in this way. Interestingly, this is 
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FIG. 1: Concurrence of the unmixed symmetric state pn (up- 
per curve), Eg. 1101 and of the pairwise mixed symmetric state 
p'n {lower curve), Eq.^O for N = 100. 



not quite correct. Rather, there are threshold values of 
n ai n = [N ± ^/N)/2, below and above which pairwise 
cancellation occurs, but between which two-state mixing 
cannot account for the loss of entanglement. Analyti- 
cally, 



C(p:J = max 



4n{N - n) 
N{N - 1) 



- 1 



(11) 



with maximal entanglement of 1/(A^ — 1) occuring at N /2 
[N even) or of 1/A'' at {N±l)/2 {N odd). Thus it is the 
least entangled states which are most resistent to pair- 
wise entanglement cancellation, as illustrated in Fig. ^ 
Further, the lack of quantum correlations in the ground 
subspace cannot be explained by incoherent mixing of the 
states within the zone {N - Vn) /2 < n < (N + ^/N) /2 
among themselves 0|. This is already apparent from 
consideration of the N = 2 and = 4 cases, where pi 
and p2 respectively have no other states in their zones 
with which to mix. But it remains true at = 6, where 
mixing of p2, ps, and p4 results in an inseparable state 
(C = 1/9). As shown in Fig. El the concurrence of the in- 
coherent mixture of the states within this zone decays as 
^ 1/N. The separability of the degenerate ground state 
therefore arises from the mixing of this inseparable core 
with the surrounding sea of pairwise-separable states. 

Whether any isotropic ferromagnet can become entan- 
gled at any temperature remains an open question. It 
seems unlikely. If the first excited state is nondegenerate, 
or nearly so, then it will change the low-temperature den- 
sity matrix p (Eq.^ by terms of order exp{— AE / kT) / N , 
where AE > is the energy difference between the 
ground and first excited state. Such changes can be seen 
to be too small to result in an inseparable state, by con- 
sidering the extreme case where the reduced density ma- 
trix for the excited state is a perfectly entangled triplet 
state. On the other hand, if the first excited state is 



FIG. 2: Concurrence of the incoherently mixed symmetric 
states in the zone (A - \/N)/2 < n < {N + ^/N) /2 as a func- 
tion of N. Within this zone pairwise entanglement cancella- 
tion fails. The curve 1/A (dashed) is shown for comparison. 



highly degenerate (as in the Heisenberg spin chain, for 
example), or if there are many closely spaced levels just 
above the first excited state, then the thermal correction 
can be of order exp(— Ai?/fcT), but in this case one ex- 
pects that incoherent mixing among the degenerate (or 
almost degenerate) excited states will result in a separa- 
ble state, as occurred when the degenerate ground states 
mixed. In either case, there appears to be no way to form 
an inseparable thermal state. 

Since the obstacle to thermal entanglement seems to 
be the highly degenerate ground state, it is natural to 
wonder if the application of an external magnetic field 
would induce entanglement by relieving the degeneracy. 
To study this question, a simple extension of a model 
introduced by Majumdar and Ghosh defined by the 
Hamiltonian 



H 



N N 

XMG — 9l ^ • Sj+i + g2 ^ Si • 
i=l 1=1 
N N 

i=l i=l 



(12) 



was studied numerically. Here _B is a homogeneous exter- 
nal magnetic field directed along the -|-z-axis, and peri- 
odic boundary conditions are assumed, so that S„j+Ar = 
Sm. The gi are coupling constants, determining the 
strength of the interaction of a spin with its ith near- 
est neighbor. When 53 = the model reduces to the 
Majumdar-Ghosh spin chain; ii §2 — Q as well, then it 
is the Heisenberg spin chain 16]. Without loss of gen- 
erality the Hamiltonian can be rescaled so that = — 1 
for a ferromagnetic chain. Concurrences were calculated 
numerically for each qubit pair in spin chains of lengths 
A^ = 4 through 8. The coupling constants (72 and 33 
were varied between —4 and 0, while the temperature 
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and magnetic field were varied between and N. In no 
case did any qubit pair exhibit any entanglement. Sim- 
ilar numerical studies of other geometries were also car- 
ried out, including open spin chains (i.e. without pe- 
riodic boundary conditions), nine spins arrranged as a 
3x3 grid with and without periodic boundary condi- 
tions, eight spins arranged as a cube, and several irreg- 
ular geometries. These calculations also failed to reveal 
any entanglement between any qubit pair. Thus it ap- 
pears reasonable to conjecture that ferromagnets cannot 
exhibit thermal entanglement at any temperature or ap- 
plied field strength. 

This paper has examined the quantum correlations 
present in isotropic ferromagnets. The ground subspace 
of an isotropic ferromagnet is spanned by the completely 
symmetric states, and it was proven that no other states 
intrude. As a result, the low-temperature two-spin re- 
duced density matrix assumes a universal, separable 
form. Pairwise entanglement cancellation cannot entirely 
account for this separability, as there exists a core of 
states which remain inseparable even when mixed inco- 
herently among themselves. Numerical evidence was pro- 
vided supporting the conjecture that isotropic ferromag- 
netic spin systems do not exhibit thermal entanglement, 
even in the presence of an external magnetic field. 

In fact, universality in the entanglement behavior of 
magnets is even more general than described here. Al- 
though anisotropic deformations of the exchange inter- 
action can modify the structure of the ground subspace 
in several ways, universal behavior does reappear at low 
temperature even in anisotropic models. This will be 
discussed in detail in a subsequent paper. 
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APPENDIX 

Theorem: The ground subspace of the A''-site con- 
nected ferromagnetic spin graph, defined by the Hamil- 
tonian H = X]i<j Jij^i'^j with Jij < 0, is spanned by 
the N + 1 completely symmetric states given in Eq. |3 
there are no 'accidental' degeneracies. 

Proof: First consider the simpler case where Jij < 0, 
that is, every spin interacts, though perhaps very weakly, 
with every other spin. The spin operator product Si • Sj 
is equivalent to ^ (Pij — ^T), where Pij is an operator 
which permutes spins i and j, and 2 is the identity op- 
erator. Each completely symmetric state is an eigen- 
state of each permutation operator separately, and so its 
energy eigenvalue is simply the sum of the coefficients 
of the individual permutation operators in the Hamilto- 
nian: Eo = (H) = i E.<, J^A'P^l " k^) = I E.<, Jv- ^ 
a state |x) is not completely symmetric, however, there 
is by definition some permutation operator which does 



not leave this state invariant. In the special case un- 
der consideration, such a permutation has a nonzero co- 
efficient in the Hamiltonian, say Jpq. Now, Vfj — I, 
hence the maximum expectation value of a permuta- 
tion operator in a normalized state is unity, so that 
Jpqixl'Ppqlx) > Jpq 7 siucc the Jij are negative numbers. 
So {x\H\x) = i E.<; JrjixlV^, - il|x) > Eo. Thus the 
completely symmetric states span the ground subspace, 
as claimed. 

Now consider the more general case, where some of the 
coupling constants Jij may equal zero, but the result- 
ing graph remains connected. Suppose once again that a 
state Ix) is not symmetric under the interchange of p and 
q. By the definition of connectedness, there is now a chain 
of spins ii, i2, . . . , in such that each of Jpi^ , Ji-^i^ , ■ ■ ■ , Ji,^q 
is nonzero. Then the (2n -I- l)th power of the 

Hamiltonian, contains a product of permutation opera- 
tors 

T^pii • • -T^in-iinT^inqT^in-lin ' ' •'^iii2^p«l (A. 13) 

which is equal to Vpq (from the right, the first n + 1 
permutations put the spin p in the position q, while the 
remaining put spin q into position p and sort the inter- 
mediate spins of the chain back into their original posi- 
tions). The above argument can now be repeated. Any 
completely symmetric state is an eigenstate of each term 
of while Ix) is not, and hence the expectation 

value of (xl-ff^"+^lx) Therefore |x) does not 

lie in the ground subspace. ■ 
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